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Abstract 

The process e’''e“ —qqgg where all hnal particles are hard is treated in a semi-analytic approach. 
Cuts on the invariant masses of the quark and the gluon pair can be applied. The total cross 
section is obtained by three numerical integrations. 


1 Introduction 

The study of W pair production is an important task of the present experiments at LEP []^. 
The W Bosons are only accessible through their decay products. Therefore, the complete process 
e+e“ —*• /1/2/3/4 must be studied. It can be calculated by Monte Carlo or semi-analytic techniques 
P. A Monte Carlo description is absolutely necessary to describe a real detector and realistic 
experimental cuts. A semi-analytic description allows only for simple cuts in the phase space but 
it is very fast and still flexible enough for hts to data. This is proved by the semi-analytic code 
GENTLE/4fan P], which is used by all four LEP collaborations to determine the W mass |l|. 

A jet due to a light quark cannot be distinguished experimentally from a jet due to a gluon. 
Hence the reaction 


e+(/ci)e (fcs) ^ q{Pi)q{P2)g{Pz)g{PA) (1) 

could give a signature similar to that of the reaction 

e+e" ^ qmqsq^- (2) 

Reaction (Q) is therefore an important incoherent background to the process (^). The process (|lD 
can be measured in more detail when the quark flavor q can be tagged. For heavy quarks, it is 
also a background to Higgs production 0. 

^email: jack@graviton.hep.physik.uni-muenchen.de 
^email: leike@graviton.hep.physik.uni-muenchen.de 
^Supported by the EC program CHRX-CT940579 




The process (|I]) was investigated many years ago p[ because a jet due to a light quark cannot 
be distinguished from a jet due to the quark accompanied by soft and/or collinear gluons. Monte 
Carlo studies of the process (|l|) with four hard final particles can be found in references p. 

The semi-analytic calculation of 4-fermion final states neglecting quark masses leads to very 
compact formulae for double and triple differential cross sections 0. Only one (two) numerical 
integrations are necessary to evaluate simple distributions (cross sections). This makes the semi- 
analytic codes fast compared to Monte Carlo programs. 

The aim of this paper is the derivation of semi-analytic formulae describing the process (P, 
which will be added to the code GENTLE/4fan in order to enhance the flexibility of this code. In 
the calculation of reaction (|I[), we must apply the same cuts as to the reaction (H). In particular, 
the invariant energies of the gluon and quark pairs in ([^) are bounded from below. This is necessary 
to keep the cross section out of the range of hadronic resonances where perturbation theory breaks 
down. Such cuts also ensure that the gluons cannot become soft. Where all four final particles 
are hard, we can limit ourselves to a tree-level calculation. 

The numerical largest corrections due to the radiation of photons from the initial state can be 
taken into account by a convolution , 

a'sys) = /■ (3) 

Js'~ S 

See reference for further details and references. 


2 Calculation 


We parameterize the phase space by 6 angles and 2 invariants. 


with 


i=l ^Pi 


2=1 


V^A(s,s„Sg) yJX{sq,ml,ml) yJX{sg,ml,ml) 


8s 


8s« 


8s„ 


dsqdsgdflodQqdQg , 


A(a, b, c) 


P\ 


+ — 2ab — 2ac — 26c, A = A(s, Sq, Sg), 

{ki + k2f, Sq = (pi P2f and Sg = (pg 

mi = m 2 = mg, m 3 = m^ = mg = 0 . 


( 4 ) 


( 5 ) 


The solid angles Vtq and Vtg are defined in the rest frames of the quark and gluon pairs, hlo is 
the solid angle of pi +'P 2 in the laboratory frame. We have dkli = d cos 6id(j)i for i = 0,q,g. The 
kinematical ranges of the integration variables are 


{mi + m2f < Sq < {^/s - m3 - m^f , {m3 + m^f < Sg < {^/s - ^/s^f, 

— 1 < cos 6i <1 , 0 < < 27r. 


( 6 ) 


The process (|^ is described by 16 Feynman diagrams, 8 with virtual photon exchange and 
8 with virtual Z-Boson exchange. 6 of the 8 diagrams are obtained by attaching the two gluons 
in all possible combinations to the quark pair, the other two diagrams contain the triple-gluon 
vertex. The two types of diagrams are presented below in Figure 1. 

We obtain distributions and cross sections by integrating over the matrix element squared, 
which is obtained by traditional trace techniques. Our semi-analytic approach does not allow 
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Figure 1: Radiation of the 2 gluons by the quarks and triple-gluon vertex 


for angular cuts to separate gluons from quarks. Therefore, we keep a non-zero quark mass to 
avoid collinear singularities. The calculation is done with the help of the symbolic manipulation 
program FORM p. 

Without transversal beam polarization, the integration over the angle 0o is trivial giving 27r. 
While the hrst 2 integrations over cfq and cfg as azimuthal angles of the qq- and ^f^f-subsystems 
in the chosen parameterization can be done relatively easily and the integration over the solid 
angle cos 6 *o is trivial because only terms proportional to 1 , cos 6^0 and cos^ arise, the last 2 
angular integrations prove to be more complicated p. This is due to logarithms with quadratic 
or biquadratic polynomials in cosOq^g in their arguments after one integration. In order to avoid 
singularities for the last 2 angles cos Oq^g - these reflect the collinear singularities of the original 
amplitude for rriq = 0 - the quark mass rriq cannot be neglected. So from these last 2 integrations 
one can be done analytically leading to dilogarithms of the type Li 2 {sg, Sg] s,mq) and the above 
mentioned logarithmic integrals as dominant contributions, while the other integration is treated 
numerically. An analytical integration over cos 6g in the massless ^f^f-subsystem with a numerical 
integration over cos Og afterwards leads to shorter expressions of the integrand than vice versa and 
is favored. 

The total cross section may then be written in the form 



^ d cos 9q ^ Cqqgg{e] q; s)ggqgg{sq, Sg, cos 9q; s, ml) 


(7) 


with numerical integrations over cos 9q, Sg and Sg. In the integration (^, we allow for cuts Sq and 
Sg on the invariant masses. After neglecting quark masses, the function Cgqgg containing couplings 
and gauge boson propagators, is 




qqggK^'i *?; ' 


Dv{s) 


2 


a^(47r)^ 3fJe 


E 

V^yj='y,Z 


1 1 
Dvi (s) Dy. (s) 


X [L,{Vi)Le{Vg) + Re{Vi)Re{Vg)] [Lg{Vi)Lq{Vq) + , 

s — My T iMy Ty. 


( 8 ) 


Lfiy) and RfiV) are the left- and right-handed couplings of the fermion / to the vector Boson V 
with the mass My and the width Ty. The preceding factor is kept for convenience with reference 
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The analytic expression for the function Qqqgg{sq, Sg, cos 9q; s,ml) is much longer than those 


obtained for the four-fermion hnal states investigated in |]^. This is partially due to the richer 


topology of the Feynman diagrams, partially due to the non-zero quark mass appearing as an 
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additional dimensional parameter and partially dne to the fact that the integration over 6q cannot 
be done analytically. As was recognized in |]^, the analytic expressions always tend to be longer 
before the integration over the last angle is done. 

Several checks are applied to ensnre the validity of onr calculation. The matrix element must 
satisfy Bose symmetry. In addition, it must vanish when the polarization vector of anyone of the 
gluons is substituted by the momentum vector of this gluon. All “artificial” poles appearing from 
partial fraction decomposition must be compensated by corresponding terms in the nominator. 
In particular, the result must be hnite in the limit A —0 and cos 6i —>■ ±1. Taylor expansions 
are necessary in different regions of the phase space to obtain results, which are numerically 
stable. All analytical integrals entering our calculation are checked numerically. In addition, our 
analytic result is checked numerically against our amplitude squared with a 7-fold integration by 


for different c.m. energies, quark masses and cuts. 

3 Results 


the integration routine VEGAS [TO. Finally, we compared our result numerically with COMPHEP |11 
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Figure 2: The total cross section & as a 
function of the c.m. energy for = 4.3 GeF 
and nic = 1.5 GeV. For > 120 GeF, the 
cross section is also shown for mb = 1.5 GeF 
and mb = 10 GeF. The cuts = y^ = 
20 GeF are applied. 


Figure 3: The differential distributions 

dajd^JTfi and da/d^JTf for ^/s = 

200 GeF, me = 1.5 GeF and mb = 4.3 GeF. 
The cuts y/s^ = y/s^ = 20 GeF are applied. 


The numerical input for our figures is a* = 0.121, a = 1/128, sin^ 9^ = 0.232, Mz = 
91.186 GeF, and Tz = 2.495 GeF. 

The total cross section for q = b,c is shown in Figure 2. To illustrate the mass dependence, 
the cross section is given for different 6-quark masses beyond the Z peak. The considered process 
has a large cross section. For light quark masses and large energies y^ > 200 GeF, it becomes 
comparable to the cross section of —>■ qq. This is due to the radiation of two co//mear gluons. 
For ^/s = 200 GeF and = 1.5 GeF, every collinear gluon generates a logarithm In^ 9.8, 

TKLq 

which numerically compensates the coupling constant Other contributions to the cross section 
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are of comparable importance for lower energies or larger qnark masses. The logarithmic mass 
terms are the reason for the increase of the total cross section with smaller qnark masses. 

Fignre 3 shows the differential distributions of the invariant energies of the quark and gluon 
pairs in the final state. The distributions are zero for y/si < 20 GeV and y/si > 180 GeV due 
to our kinematical cuts. The distributions for peak for small Sg reflecting the infrared peak 

due to the radiation of soft gluons. The distributions for are largest for values of Sg near the 

upper limit (^i — ^/s^y again due to the radiation of soft and collinear gluons. 

To summarize, we performed a semi-analytic calculation of the process —>■ qqgg where 

all final particles are assumed to be hard. The process has a large cross section. It is a large 
incoherent background to the four-fermion processes currently measured at LEP. 
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